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1 . FOREWORD 


The Solar Energy System Economic Evaluation - F'-^al Report has been 
developed by the George C. Marshall Space Flight Center as a part of 
the Solar Heating and Cool 'ng Development Program funded by the 
Department of Energy. The analysis contained in this document de- 
scribes the economic performance of an Operational Test Site (OTS). 

The objective of the analysis is to report the long-term economic per- 
formance of the system at its installation site and to extrapolate to 
four additional locations which have been selected to demonstrate the 
viability of the design over .'4 broad range of environmental and economic 
conditions. 

The contents of this document are divided into the following topics: 

• System Description 

t Study Approach 

• Economic Analysis and System Optimization 

t Results of Analysis: Technical and Economic 

f Economic Uncertainty Analysis 

• Summary and Conclusions 

The data used for the economic analysis have been generated through eval- 
uation of the Operational Test Site described in this document. The data 
that have been collected, processed, and maintained under the OTS Develop- 
ment Program provide the resource from which inputs to the simulation 
programs used to perform technical and economic analysis are extracted. 

The Final Report document, in conjunction with the Seasonal Report for 
each Operational Test Site in the Development Program, culminates the 
technical activities which began with site selection and instrumentation 
system design in April, 1976. The Seasonal Report emphasizes the technical 
analysis of solar systems performance. It compares actual performance 
with predicted performance derived through simulation methods where 
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actual weather and loads defined the inputs. The simulation used for 
final report analysis is based on the technical results of the seasonal 
report simulation, with the exception that long-term weather, and derived 
loads are used as inputs instead of measured weather and loads. This 
causes the expected value c'f solar system performance in the Seasonal 
and Final Reports to differ. In addition, localized and standard eco- 
nomic parameters are used for economic analysis in the final report 
evaluation. The details of the simulation program are described in 
Reference [1].* Other documents specifically related to the solar 
energy system analyzed in this report are [2], [3] and [7]. 


■^Numbers in brackets designate references found in Section 8. 
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2. SYSTEM DESCRIPTION 


The IBM-System 2 Solar Energy System is located in a single-family 
dwelling at the Veterans Administration Center in Togus, Maine. The 
system is designed to preheat approximately 56 percent of the domestic 
hot water.' Silicon fluid is circulated through a 105 square-foot 
Libbey-Owens-Ford (L.O.F) flat-plate collector array, and a double wall 
heat exchanger. The collector array faces 15 degrees west of due so'th, 
at an angle of 45 degrees from the horizontal. Water from the 120- 
gallon preheat tank is circulated through the other side of the heat 
exchanger. The preheat tank services a standard electric 40-gallon 
domestic hot water CDHW) heater, which adds the necessary auxiliary 
energy. Figure 2-1 is a schematic of the system. Sensor designations 
are in accordance with NBS-IR-76-1137 [4]. The measurement symbol 
prefixes: W, T, EP, and I represent respectively: flow rate, tempera- 
ture, electric power, and insolation. Figure 2-2 shows two pictorial 
views of the installation. This system has only one mode of operation. 
The system turns on when the collector outlet temperature becomes approx- 
imately 25°F above the temperature of water in the bottom of the preheat 
tank, and turns off when this differential temperature falls below 8°F. 
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Figure 2-1. IBM System 2 Schematic 




OWC.TNM 


CoUector Array During Installation 


DHW Heater and Preheat Tank in Basement 
Figure 2-2. IBM System 2 Pictona! 





3. STUDY APPROACH 


3.1 Introduction 

The Final Report is an economic evaluation of the solar energy system 
(based on life cycle costs versus energy savings) for five cities which 
are considered to be representative of a broad range of environmental 
and ecomonic conditions in the United States. Life cycle costs provide 
a measure of the total costs of owning and operating a system over the 
life of the system rather than focusing solely on the initial cost of the 
system. The life cycle costs used in this evaluation consider hardware, 
installation, maintenance, and operating costs for the solar-unique 
components of the total system. Energy savings result from replacement 
of conventional forms of energy by solar energy after the costs of pro- 
ducii’g the solar energy ai*s deducted. The total system operates in a 
scenario that comprises long-term average environmental conditions, 
loads, fuel costs and other economic factors that are applicable in 
each of five cities. 

The five cities include four standard analysis sites which were selected 
according to the criteria listed below and the site where the system 
was, in fact, installed and operated. The selection cr'-teria were based 
on: 


• Availability of long-term weather data 

t Heating degree days (load related factor) 

• Cold water supply temperature (load related factor) 

• Solar insolation 

• Utility rates 

• Market potential 

• Type of solar system 

To achieve the range of environmental and economic parameters desired, 
the four locations listed below, plus the actual installation location, 
were used. A solar energy system buyer may evaluate his own local 
environmental and economic conditions relative to those considered in 
this Final Report by comparing the insolation available, the heating 
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load (applicable only to space heating systems), and the utility rates 
against the results reported in Section 5. 

Albuquerque, New Mexico 

1828 Btu/Ft^/day average insolation* 

Medium heating load (4292 HDD) 

High utility rates (>0.06 $/kWh)** 

Fort Worth, Texas 
2 

1475 Btu/Ft /day average insolation* 

Light heating load (2382 HDD) 

Medium utility rates (0,04-0.06 $/kWh)** 

Madison, Wisconsin 
2 

1191 Btu/Ft /day average insolation* 

High heating load (7730 HDD) 

Medium utility rates (0.04-0.06 $/kWh)** 

Washington, DC 

2 

1208 Btu/Ft /day average insolation* 

Medium heating load (5010 HDD) 

High utility rates (>0.06 $/kWh)** 

Actual Site (Tonus, Maine) 

n 

1228 Btu/FtVday average winter insolation* 

High heating load (7946 HDD) 

Medium utility rates (0.04-0.06 $/kWh)** 

The parameters that define the system design were derived from the actual 
operating conditions of the system at the installation site. Solar energy 
system design may be economically optimized for the site at which the 


*Insolation values are average daily long-term values on a horizontal 
surface. 


**Ut11ity rates are effective yearly average values based on 1(XW vse 
for schedules in effect for January 1980. See Appendix D. 
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system Is Installed. The fundamental abjactlve In optimizing the design 
of a solar energy system on an economic basis is to minimize cost by 
allocating the required amount of energy between the solar and conven- 
tional portions of the system. To attain this objective, each unit of 
energy should be produced by the portion of the total system which 
generates the lowest incremental cost in producing that additional unit 
of enirgy. This e accomplished in the final report analysis by determin- 
ing the optimal solar energy system size (collector area or equivalents 
ly, solar fraction). 

In the Operational Test Site (OTS) Development Program there are many 
solar energy systems designed by many different contractors. Some of 
the designs were installed in new buildings and some were retrofitted to 
existing buildings. Consequently, there are a variety of factors which 
contributed to the design of a system at a given site. In some cases 
the objective of optimizing the design according to the previously 
stated criterion could not be met. A method of evaluation which estab- 
lishes a common basis for evaluation of all these systems was required. 
The method selected is to optimize the collector size through the 
f»Chart [1] design procedure. F-Chart is a design program developed 
by the University of Wisconsin for solar heating and/or domestic hot 
water systems. The program uses a set of design charts (developed 
by detailed simulations) which estimate the thermal performance of 
a solar system based on collector characteristics, storage, energy 
demands, and regional long-term weather data. Using the results 
of thermal analysis, an iterative procedure is implemented to select a 
collector area which minimizes the life cycle costs. Once the optimal 
collector size has been determined, the resulting thermal and economic 
performance can be obtained. 

The resolution of two inter-reUted problems was required in order to 
adapt f-Chart to the evaluation developed in the Final Report. The 
first was how to use the data and experience gained from the actual 
operation of the solar energy system; the second was what procedure to 
follow in view of the fact that all solar energy systems to be analyzed 
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do not have optimal collector area sizing. To resolve the first problem, 
the characteristics of design and operation of the existing solar energy 
system were used to develop the input parameters for f-Chart. This 
procedure, detailed in Appendix A, involved the normalization of collec- 
tor flow rates and storage capacity to collector area. Collector charac- 
teristics developed from field data through a collector analysis program 
were substituted for the theoretical single panel parameters furnished 
by the collector manufacturer. To resolve the problem of different 
collector areas, an optimal collector area was derived for the site. 

The final adaption of f-Chart includes the Inputs derived from opera- 
tional data and optimal collector area. 

As the system application at each of the five analysis sites is studied, 
the loads are iteratively redefined, the site peculiar- parameters are 
changed as described in Appendix A, and a new optimal collector area is 
computed. The economic factors are the result of the f-Chart analysis 
with these inputs. 
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3.2 Groundrules and Assumptions 


The cost differential between solar and the conventional system Is 
significant to the economic evaluation In the Final Report. Cost Items 
which were equal for both alternatives do not contribute to the differential 
cost. The cost of the conventional system was assumed to be Identical 
with or without the solar alternative. Although a conventional system 
Is usually selected according to the availability and cost of energy In 
a particular geographic region, this alternative Is not permitted In the 
final report analysis because an existing system is being evaluated. 

Savings which might be realized by comparing solar against an auxiliary 
other than the design option were not evaluated. The system configura- 
tion, Including the conventional auxiliary. Is the same for all five 
analysis sites. 

The cost of the solar unique hardware is based on mass production es..-- 
mates. The total Incremental costs for acquisition of a solar alterna- 
tive are tha sum of a cost proportional to collector area and a cost 
Independent of collector area. For economic evaluation, life cycle 
costs (I.e., costs of acquiring, operating and maintaining the solar 
systems) were forecast on an annual basis over the design lifetime of 
the system, then discounted to an equivalent single constant dollar 
(1980) value as described In Section 4. 

Fuel costs are calculated at current (1980) local values for each of the 
five analysis sites. Other economic parameters are standardized by 
referencing current national economic conditions. Maintenance, Insur- 
ance, depreciation, system life, salvage values (for commercial systems) 
are determined from best experience. Tax credits allowed by the Federal 
Government for the solar energy systems are credited against the acqui- 
sition cost. A combined state and federal income tax rate of 30 percent 
Is assumed for estimating tax savings resulting from the Interest paid 
In financing a solar system. Property taxes arising from the Increased 
value of property with an Installed solar system are neglected due to 
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the current trend in many states to forego these taxes to prevent them 
from being a hindrance to solar energy usage. 

The primary measure of cost effectiveness of the solar energy system 
in the Final Report is: 

• Life Cycle Cumulative Savings (LCCS) - The present value of 
the cumulative energy savings (in dollars) that result from 
operation of the solar system alternative instead of the 
conventional backup. 

Two secondary measures that depend on the life cycle cumulative savings are 

• Year of Positive Savings - Year in which solar system first 
becomes profitable; i.e., the annual conventional fuel bill 
without solar exceeds the sum of the annual fuel bill with 
solar and the annual cost for the solar system. 

f Year of Payback - Year in which the compounded net savings 
equals the initial cost for the solar system. Net savings 
were computed with respect to the fuel cost of the conven- 
tional system. 
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4. ECONOMIC ANALYSIS 

4.1 Factors in Life Cycle Costs and Savings 


The economic calculations of this study are performed in the f-Chart 
program and are based on comparisons uf l^fe cycle costs of conventional 
energy systems with those of solar erer;* • systems. The life cycle sav- 
ings of a solar energy system over a conventional energy system can be 
expressed as the difference between the total fuel savings that result 
from operation of the solar energy system and the increased costs that 
result from the Investment in, the operation of, and the maintenance of the 
solar energy system. The savings can be expressed by the relationship: 

LCCS = (Cp/n)LF - (C^A + C^) (1) 

where LCCS = Life cycle cost savings of the solar 

energy system ($) in terms of present worth 

p.j = Factor relating life cycle fuel cost savings 
to first year fuel cost savings 

Cp/n = Fuel cost per unit divided by conventional 
heating unit efficiency 

L = Total load on system computed from long- 
term average conditions (Btu) 

F = Solar fraction 

?2 ~ Factor relating life cycle investment 
operation and maintenance expenditures 
to the Initial investment 

= Solar energy system costs dependent 
on the collector area ($/Ft^) 

A = Collector area (Ft^) 

Cp = Solar energy system costs that are independent 
of collector area. ($) 
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It Is assumed that the costs of components which are common to both 
conventional and solar heating systems (e.g. the furnace, ductwork, 
blowers, thermostat, etc.), and the maintenance costs of this equipment, were 
identical. Consequently, all references to solar energy system costs 
refer to the cost Increment above the common costs. 

The multiplying factors, and Pg, facilitate the use of life cycle 
cost methods in a compact form. Any cost which is proportional to either 
the first year fuel cost or the initial investment can be included. These 
factors allow for variation of annual expenses with inflation and they 
reflect the time value* of money by discounting future expenses to present 
dollar values. 

To illustrate the evaluation of P.| and P., consider a simple economic 
situation in which the only significant cobts are fuel and system equip- 
ment costs. The fuel cost is assumed to escalate at a constant annual 
rate, and the owner pays cash for the system. Here, P.| accounts for fuel 
escalation and the discounting of future payments. The factor P^ accounts 
for investment related expenses which, in this case, consist only of the 
Investment which is already expressed in current dollars. The factors P^ 
and Pg are then 

P^ = PWF(N, e, d) 


where N ■ Period of economic analysis (yrs) 
e ■ Escalation rate of fuel price 
d ■ Annual discount rate 

*D1scounting refers to the fact that an expense that is anticipated to be 
$1000 in 10 years is equivalent to an Investment today of $463 at a discount 
rate of 8%. 
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The function PWF(N, e, d) is the present worth factor that accounts for 
inflating payments in discounted money. 


PWF(N. e, d) = -g-l-i 



( 3 ) 


When multiplied by a first period cost (which is inflated at a rate, e, and 
discounted at a rate, d, over N years), the resulting value is the present 
worth life cycle cost. 


In the more complex analysis the expenditures incurred by the additional 
capital investment cause and to take the following form: 


P^ = (1 - Ct) PWF(N, e, d) 

^2 *^21 ^22 " *^23 *^24 ''' ^ 2b “ ^26 " ^27 

where Pg^ = Factor representing the down payment 

Pg2 = Factor representing the life cycle cost 
of the mortgage principal and interest 

P23 = Factor representing income tax deductions 
for interest payment 

= Factor representing miscellaneous costs 
(maintenance, insurance, etc) 

P25 = Factor representing net property tax costs 


( 4 ) 

( 5 ) 


P2g = Factor representing straight line depreciation 
tax deduction for commercial installations 


P27 = Factor representing salvage (commercial installation) 
or resale value (residential installation). 
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The factors 


where 


P 21 through P 27 are defined as follows: 


21 “ 

( 6 ) 

'22 =* (1 “ D) PWF(N, 0, d)/PWF(N, 0, 1) 

(7) 

'23 (1 - D) t)pWF(N, i, e) [i - 1/PWF(N, 0, i)] 

. +PWF(N, 0, d)/PWF(N, 0, i)[ 

( 8 ) 

'24 = (1 - Ct) M PWF(N, g, d) 

(9) 

'25 = t (1 - t) V PWF(N, g, d) 

( 10 ) 

'26 = (Ct/N) PWF(N, 0, d) 

01 ) 

27 = G/(l + df 

( 12 ) 


D = Ratio of down payment to the initial investment 

N = Period of analysis (Note that the period of analysis, 
the term of the loan, the depreciation lifetime, and 
the years over which the depreciation deductions con- 
tribute to the analysis are arbitrarily set equal in 
this report). 

d = Discount rate (after tax return on the best 
alternative investment) 

i = Annual mortgage interest rate 

t = Effective income tax rate 

C = Commercial or non-commercial flag (1 or 0 
respectively) 
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M = Ratio of first year miscellaneous costs to 
Initial Investment 

g = General Inflation rate 

t = Property tax rate based on assessed value 

V = Ratio of assessed value in first year to initial 
investment 

G = Ratio of salvage or resale value to initial 
investment 

For a given location, heating load, and economic situation, it is possible 
to optimize the system design variables to yield the maximum life cycle 
savings. The main solar energy system design variable is the collector 
area. The effect of collector area on the life cycle savings is illus- 
trated in Figure 4-1 for the four sets of economic conditions. Curve A 
corresponds to an economic scenario in which solar energy cannot compete 
with the conventional system. Curve B exhibits a non-zero optimum area, 
but the conventional system is still the most economical. Curve C 
corresponds to the critical condition where solar energy can just compete 
with the conventional system. Curve D corresponds to an economic scenario 
in which the solar energy system is the most economical. 

Each curve of Figure 4-1 begins with a negative savings for zero collec- 
tor area. The magnitude of this loss is C^, and reflects the presence 
of solar energy system fixed costs in the absence of any fuel savings. 

As the collector area increases Curves B, C, and D show increased sav- 
ings until reaching a maximum at some optimum collector area. As the 
collector area is further increased, the fuel savings continue to increase, 
but the excessive system cost forces the life cycle savings of the sys- 
tem to decrease. The collector areas at each of the five analysis sites 
listed in this report have been optimized by the f-Chart program analysis 
technique for the long-term average weather conditions and the economic 
conditions at that site. 
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4 . 2 Federal Tax Credits for Solar Energy Systems 

The Federal Government has provided tax incentives that are applicable 
to solar energy systems.* This credit is 30 percent of the first $2000 
plus 20% of the next $8000 spent on solar equipment, or a maximum credit 
of $2200. The credit is applied in this analysis by reducing both the 
collector area dependent cost and the cost independent of the collector 
area, or constant solar cost, by an effective credit factor based on the 
total cost of the system, 

As an example of the tax credit computation, assume the collector area 
dependent cost is $30/ft based on 100 ft and the constant solar cost 
is $900 for a total price of $3900. The effective credit factor is: 

2000 X 0.30 + (3900 - 200 0) x 0.2 0 _ « o 
3900 

Therefore the adjusted costs used as f-Chart inputs are: 

Collector area dependent cost 
C^, = $30 X (1 - 0.2) = $24/ft^ 

Constant solar cost 

Cg, = $900 X (1 - 0.2) = $720 

The f-Chart econonomic analysis is modified by using these adjusted 
costs to reflect tax credit effects. Including tax credit in area 
optimization is an iterative process since the credit is affected by the 
system size and vice versa. Optimal collector area was modified in this 
analysis, as were the f-Chart economic parameters, by use of the tax 
credit. Items 23 and 24 in Table 5.1-2 reflect the solar costs before 
application of tax credits in terms of collector area dependent cost and 
constant solar cost. Initial system costs before and after tax credit 
inclusion are shown in Table 5.2-1 for each site based on optimal col- 
lector area. 

* The tax credit has been revised after 1979 to 40 percent of the first 
$10,000 for a maximum credit of $4,000. The new effective credit factor 
as given in the example above is 0.4 for systems o .ting less than $10,000 
and the ratio of the maximum credit to the total s> ^tem cost for systems 
costing more than $10,000. 


I 
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5. RESULTS OF ANALYSIS 


5,1 Technical Results 


For each of the five analysis sitns ah optimal solar system based on the 
con.iguration of the actual 1nsta‘*1ation was determined by using the f- 
Chart design procedure, The environmental parameters and the loads used 
in this procedure for each of the five sites are shown in Table 5.1-1, 

In applying the design procedure, a process that iterates on the col- 
lector area was used. Figures 5.1-1 (a)-(e) show the results of that 
design procedure in terms of the expected solar fraction versus the 
collector area for each site. The expected solar fraction is the ratio 
of the expected solar energy used toward satisfying the load to the 
total load. The graphs in Figures 5.1-1 (a)-(e) show that as the collectcr 
area increases, the expected solar fraction ‘increases asymptotically. 
However, the economically optimal collector area was selected to maximize 
the economic benefits of the solar energy system, not the expected solar 
fraction. The optimal collector area is shown by the dotted line for 
each site. Increasing the collector area beyond the optimal value 
forces a diminishing return on the investment for the system. The 
expected solar fraction for the optimal collector area is shown in the 
last column in Table 5.1-1. 

The resulting thermal performance, once the optimal size system is se- 
lected, is shown in the graphs of Figures 5.1-2(a)-(e) for each analysis 
site. The incident solar energy was derived from long-term average inso- 
lation at the site. The total load was computed based on design parameters 
of the actual system as installed, modified by environmental conditions 
at each site. The load calculations are detailed in Appendix A. The 
useful solar energy is the product of the system solar fraction and the 
total load and shows on a month by month basis the portion of the total 
load that is expected to be supplied by solar energy. The shaded por- 
tion between the total load curve and the curve of useful solar energy 
must be supplied by conventional energy. 
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TABLE 5.1-1 
SUMMARY TABLE 


SOLAR SYSTEM LOAD FACTORS AND ENVIRONMENTAL PARAMETERS 



' TOTAL ANNUAL LOAD 

— T 



1 



i (MILLION BTU) 

(1 

1 ENVIRONMENTAL PARAMETERS - LONG-TERM 

& 

t 

1 

1 

1 

1 

j 



HOT 

1 INSOLATION 

HEATING 

SUPPLY WATER 

i 

EXPECTED SOLAR 

SITE 

! 

WATER 

^ BTU/FTVDAY 

DEGREE DAYS 

TEMP (°F) 

! 

_L 

FRACTION* j 


* 


TOGUS 

22.8 

1228 

7946 

52 

! 75.4 

fl 

ALBUQUERQUE 

17.2 

1828 

4292 

73 

\ 93.1 

FORT WORTH 

19.36 

1475 

2382 

65 

a 

« 73.6 

!! 

' MADISON 

22.1 

1191 

7730 

54 

j 

: 51.7 

WASHINGTON 

20.7 

1208 

5010 

60 

77.8 

’i 


*For optimal collector area 
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Table ^/.Vs 

0 „ D 0 

0.v(>9 BTU/H-r I 

0 * 0 

. 2*00 

fablT5.1-3' 

Table 5.1-3 

9„4':) BUJ/L "FIP 

Table 5.1-3 

Table 5.1-3 

AV. 00 (vAI./DAY 

1.''6,>10 F 

Table 5.1-3 

4- * 00 

1 ,^0A 
[ * 00 
2.00 

0*0 %/YR 
*;iO*00 YEaR,Y 
12.* 70 III /FT*'/ COLl 
1413*00 !|i 
2.0*00 X, 

13*00% 

20*00 YEAR,? 

6..,!.) 0 2V 

0*50 % 

i 0*00 % 

Table 5.1-3 

fftOO 

HlTtj?)’* ;t ’ 

See Note 1 

1*00 

12*l:i'0 % 

1 * 00 

;T£.*00 % 

0'jr~ 

6*00 % 

2*00 
0 * 0 % 

2 * 00 




NOTE: 1. Since the backup for the solar system is assumed to be the same type 
of system as would conventionally be used without a solar system, 
backup fuel costs and conventional costs per million Btu are equal. 
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ALBUCUERQUE, NEW MEXICO 
OPTIMAL COLLECTOR AREA - 168 FT2 



COLLECTOR AREA (FT^) 


Figure 5.1-1 (bj Solar Fraction vs Collector Area for Albuquerque, New Mexico 
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Figur$ 5.1 -1 (c) Solar Fraction vs Collector Area for Fort Worth, Texas 
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WASHINGTON. D. C. 

OPTIMAL COLLECTOR AREA » 252 FT^ 



Figura 5.1-1 {a) Solar Fraction w Collector Area far Washington, D. C. 
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TOGUS, MAINE 

OPTIMAL COLLECTOR AREA - 273 FT* 









\ 







/ 



■ 



\ 






( 






\ 

V 




j 


3IDENT 

ERGY 

SOLAR 





\ 




/ 

E^ 





1 




/ 









m 


TA 

1 

■ 


m 

1 

m 

1 

n 


n 

1 

n 

1 

1 

■ 

1 

■ 

■ 

■ 

■ 

■ 


a 















ED ARE 
JXILIAR 

^ TO Bl 
Y ENER 

E SUPPL 
GY 

“USEF 

lED 





1 

a 


* 

JL SOL/ 

•OTAL 1- 
\R ENEf 

n 


■ 

B 

B 



MONTH - 1979 


oRiGiNAi: PAGE rs 
P.E POOR, .QUALITY 


figure 5.7-2 (e) Thermal Performance of Solar Energy System with Optimized Collector Area for Togus, Maine 
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ALBUQUERQUE, NEW MEXICQ 



fijMr* S.1-2 (b) Thtrmal Pmformnca of Solar Energy System with Optimized Collector Area for 
Albuquerque, New Mexico 
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FORT WORTH, TEXAS 



Figure 5.1-2 (c) Thermal Performance of Solar Energy System with Optimized Collector Area for 
Fort Worth, Texas 
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WASHINGTON, D. C. 

OPTIMAL COLLECTOR AREA - 262 FT* 
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Figure 5.1-2 (e) Thermal Performance of Solar Energy System with Optimized Collector Area for 
l/^shington, 0. C. 
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The optimum collector for the five sites varied from 147 square feet at 
Madison to 273 square feet at Togus. The Albuquerque site with 168 square 
feet should provide more useful solar energy throughout the year than the 
other sites. Due to relatively low fuel cost and relatively poor insolation 
at Madison, the optimum collector area was the lowest of the five sites. 
Togus, with high fuel cost and average incolation, could justify the largest 
collector area. 

The technical parameters that describe the solar energy system are 
listed in Table 5.1-2 as Items 1 through 21. These parameters are 
described in detail in Appendix A. Their values are listed by site 
in Table 5.1-3. The remaining technical parameters are assigned values 
which are constant for all sites. 

The economic parameters for the solar energy system are listed in Table 
5.1-2 as Items 22 through 43, and are also described in Appendix A with 
the source for the assigned value designated. 

The following items are a function of the analysis site. 

• Collector area 

• Collector slope 

• Azimuth angle 

• Effective building DA (applicable to solar space heating systems) 

• Water main temperature 

• Present cost of solar backup fuel 

• Present cost of conventional fuel 

These are listed by site in Table 5.1-3. 
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SOLAR SYSTEM TECHNICAL PARAMETERS FOR F-CHART PROGRAM 



*An effective rate is computed for each location based on 1000 kWh use. This effective rate includes all 
charges specified in the rate schedules in Appendix D. 


5.2 Economic Results 


An essential factor- in maximizing the life cycle savings of a solar 
energy system, or conversely, of minimizing life cycle costs Is the 
economic optimization of the collector area based on equipment and fuel 
(conventional energy) costs and the capability of the solar system to 
replace significant quantities of conventional energy with solar energy. 

The replacement capability is directly dependent on the environmental 
conditions at the installation site, i.e. available solar energy. 

The graphs of Figures 5.2-1 (a)-(e) show the relationship of the factors 
comprising life cycle costs - equipment costs and fuel costs - as a 
function of collector area. Both costs are presented in terms of present 
value, i.e. baselined to today's dollars. It can be readily seen that 
as collector area increases, solar equipment costs increase proportion- 
ately. Also, as collector area increases the fuel costs decrease, 
although not as a straight line function. At some given collector area, 
the sum of these two costs is a minimum, as shown by the life cycle cost 
(LCC) curve. This minimum defines the optimal collector area for the 
given installation site. 

The solar equipment costs discussed in the preceding paragraphs include the 
principal and interest paid on an assumed 13,5 percent, 20 year mortgage, the 
income tax deduction for interest for an owner in the 30 percent bracket 
and the insurance and maintenance costs estimated at 0.5 percent of the 
initial costs. The fuel cost is that which is required by the conven- 
tional backup system and includes the effects of the f-Chart solar 
system model . 

The life cycle costs are not to be confused with life cycle savings. 

Life cycle savings is the difference between the life cycle cost of 
fuel for a conventional system and the life cycle cost of owning, 
operating and maintaining a solar energy system. 
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A summary of the costs and savings for the conventional system and the 
solar energy system is shown in Table 5.2-1 in terms of today's dollars 
expended over the analysis period. It should be recalled that the 
equipment costs shown do not include the cost of the conventional system 
since this system must be provided with or without the solar energy 
system, The equipment costs include only the additional hardware that 
must be provided for the solar energy system. This includes the 
following: 


• Collectors and mounting hardware 

• Piping and duct work (including valves and dampers) 

• Heat exchanger(s) 

• Storage un1t(s) 

• Control system 

The best estimates of equipment costs for solar energy systems indicate 
that costs fall into two categories » (1) costs dependent on collector 
area and, (2) costs independent of collector area (constant), This is 
the case, especially for residential systems, because regardless of the 
exact collector area used, certain items of equipment must be provided 
and the costs of hardware and labor for installation seem to be rela- 
tively constant. However, the cost of collectors, and certain incre- 
mental costs, are dependent on the size of the collectors used. These 
costs are shown in Table 5.2-1 for each of the five analysis sites and 
the total cost for the system is the sum of the constant cost and the 
area dependent cost multiplied by the collector area. 

The initial cost of the system in this analysis should be adjusted for 
the federal tax credit (and any other tax credit allowed by the state or 
local governments) by the methods discussed in Section 4.2. These ad- 
justed costs are shown in parentheses under "Initial Cost of System" in 
Table 5.2-1 and are used in computing the "Present Worth of Total Solar 
Costs." 
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TOGUS, MAINE 

OPTIMAL COLLECTOR AREA >• 273 FT^ 
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fl/m 5J-1 Ot) Optimiation of Collector Area for Albuquerque, New Mexico 
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Fiaun 5.2-1 fc) Optimization of CoHector Area for Fort Worth, Texas 
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MADISON, WISCONSIN 



Figun SJ-1 (d) Optimization of Collector Area for Madison. Wisconsin 
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WASHINGTON, D.C. 



Figure 5.2-1 (e) Optimization of Collector Area far Wadtington, D.C. 
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Some conventional energy must be expended with or without the solar 
energy system because, in most cases, the solar energy system will 
replace only a portion of the total energy required to support the load. 
Savings are possible with the solar energy system only because the total 
costs with the solar energy system are less than the costs of conventional 
energy. Consequently, the fuel costs over the analysis period (20 years) 
are shown in Table 5.2-1 with and without the solar energy system. 

It is assumed in this analysis that the solar energy system would be financed 
through a 20 year loan at an interest rate of 13.5 percent. Property 
taxes are assumed to be zero, but this may not be universally true. 

Insurance on the value of the solar energy system and maintenance costs 
are assumed to be 0.5 percent per year of the initial costs. Since 
interest paid on a loan is tax deductible for federal taxes, and in most 
cases for state taxes, at different rates according to the income tax 
bracket of the borrower, a 30 percent combined federal -state tax bracket 
was assumed. The value of all these costs based on the assumptions of 
this analysis is shown as the "Present Worth of Other Solar Costs" in 
Table 5.2-1. Combined with the costs for fuel with the solar energy system, 
the value is the "Present Worth of Total Solar Costs." 

Since only incremental equipment and associated costs are included in 
the analysis, the present worth of total costs for the conventional 
system without solar are simply the cost of fuel without solar. Then 
the "Present Worth of Cumulative Savings" is the difference between the 
"Present Worth of Total Costs Without Solar" and the "Present Worth of 
the Total Costs With Solar". These values for each of the five analysis 
sites are listed in Table 5.2-1. 

Finally, two economic performance parameters called "Year of Positive 
Savings" and the "Year of Payback" are shown in Table 5c2-l. As previ- 
ously discussed the year of positive savings is the year after purchase 
in which the solar system first becomes profitable, i.e., the annual 
conventional fuel bill without solar exceeds the sum of the annual fuel bill 


SUMMARY TABLE 
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TABLE 5.2-1 

COSTS AND SAVINGS OVER 20 YEAR ANALYSIS PERIOD IN DOLLARS (1980) 


OPTIMAL COLLECTOR AREA 


I 

SITE 

INITIAL COST OF SYSTEM^ 

PRESENT WORTH 
OF FUEL COSTS 

PRESENT PRESENT 

WORTH WORTH 

OF OF 

OTHER TOTAL 

SOLAR i SOLAR 

COSTS j COSTS 

PRESENT 
WORTH 
OF TOTAL 
COSTS W/0 
SOLAR 

J 

i 

PRESENT i 
WORTH i 

OF YEAR OF 

CUMULATIVE POSITIVE 
SAVINGS SAVINGS 

YEAR OF 
PAYBACK 

CONSTANT 

AREA 

DEPENDENT 

TOTAL 

WITH W/0 

SOLAR SOLAR 

TOGUS 
273 FT^ 

■ 

4641 
(3481 ) 

6525 

(4894) 

1 

2298 9371 

t 

5703 ! 8001 

! 

: 

1 

9371 

Ir 

1369 6 

- 
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ALBUQUERQUE 
168 FT^ 

■ 

2856 

(2142) 

1 , j j j , 

4740 ; 639 9318 : 4143 : 4782 9318 * 4536 1 12 

t , g 

(3555) : ; ' ’ 

FORT WORTH 
147 FT^ 

m 

2499 

(1874) 


MADISON 
147 FT^ 

m 

2499 

(1874) 

4383 i 

i 

(3287) 

i 1 

3119 7183 i 3833 ' 6952 j 7183 

' i 

; ‘ i 

231 8 1 >20 

i ( 

; t 

WASHINGTON 
252 FT^ 

m 

4284 

(3213) 

61 68 
(4626) 

2412 10886 ; 5390 * 7802 * 10886 ! 3038 4 

: 1 1 i 

^ ^ ^ ^ : 

3 

1 15 

L 


NOTES: 


I 1, Values in parentheses are approximated for the Federal tax credit by the method detailed in Section 4.2. 
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with solar and the annual costs for the solar system. The year of pay- 
back is the year after purchase when the compounded net savings equals 
the initial investment for the solar system. The factors that determine 
years until positive savings are shown in Figures 5.2-2(a)-(e) for each 
analysis site. The factors that determine the years until payback are 
shown in Figures 5.2-3(a)-(e) for each analysis site. 

A study of the Life Cost Cycle Cost curves in Figures 5.2-l(a)-(e) shows 
that an optimal collector area exists for all 5 sites. Further economic 
analysis is required to demonstrate the economic feasibility of install- 
ing a System 2 type solar energy system at each of the analysis sites. 

From Figures 5.2-2(a)-(e) , the annual expense for heating domestic hot 
water is shown to be less with solar. The additional information from 
Figures 5.2-3(a)-(e) is required to show that a system of this type is 
not economically feasible for Madison; although, the system would pro- 
vide a positive savings within the 20 year life at the other sites. 


TOGUS, MAINE 



Fi§un 5.2-2 fa) Annua) Bxpanm for Solar System and Conventional System 
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Fifun SJ2-2 (i) Annual Expanus for Solar Systam and Conventional System 
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ANNUAL COST IN DOLLARS 


FORT WORTH, TEXAS 

OFTIMAL COLLECTOR AREA « 147 FT^ 



Figun 5.2-2 (c) Annut! Expenm for Solar System and Conventional System 


MADISON, WISCONSIN 



Fi/un SJ-2 (d) Annual Expanses for Solar System and Conventional System et Madison, Wisconsin 
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WASHINGTON, D.C. 

OTTIMAL COLLECTOR AREA - »2 FT* 



Figuf* 5.M (t) AnrHM/ Ex/tmms for Solar Systam and Conventional System 
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TOGUS, MARINE 



Fifun 5.2-3 (a) Payback for Solar Energy System for Tagus, Mair~ 
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ALBUQUERQUE, NEW MEXICO 
OPTIMAL COLLECTOR AREA - 1M FT^ 



YEARS 


Figurt 5JF3 (b) Ptybeck for Solar Enargy Systam for A/biiquerquo, N'aw Mexico 
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Fifun SJ-3 (c) Ftybtck for Solar Enargy System for Fort Worth, Texet 
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MADISON, WISCONSIN 



Fi ftn 5J-3 (d) PtybMck for So/ar Bnergy System for Madison, Wisconsin 
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6. ECONOMIC UNCERTAINTY ANALYSIS 

The economic evaluation methods presented in this report are based on 
the assumption that reliable values for economic variables can be as- 
signed, However, there is an inherent uncertainty in predicting future 
expenses and benefits which is magnified by international economic 
unstability. As a consequence, the results of both the life cycle cost 
analysis and the optimization procedures must be accepted with dis- 
cretion and the effect of uncertainties must be evaluated. 


For a given set of conditions, the change in the present worth of life 
cycle cumulative savings (Table 5,2-1), aLCCS, resulting from a change 
in a particular variable. Ax-, can be approximated by the following: 

w 


_ 8LCCS 

ALCCS *“ Ax • 

axj 3 


(13) 


The expression for 8LCCS/ax^ can be obtained by direct differentiation 
of the life cycle savings equation. The life cycle cost model of 
Equations (1), (4) and (6)-(12) will be used for this analysis. The 
derivatives of these equations for each variable are given in Appendix B. 

To illustrate the use of these relationships, Uncertainty Analysis 
Tables 6-1 through 6-5 were made up for the installation site. The 
tables give the change in solar system life cycle cumulative savings, ALCCS, 
caused by a 10 percent relative increase in each of the variables. 


Table 6-1 for Togus, Maine, shows, for example, that a 10 percent in- 
crease in the discount rate from 8.5 to 9.4 percent yields a decrease in 
the value of P-j of approximately 2.434, giving a modified value of = 
24.136 (9.2 percent decrease). The value 6f decreases by 0.069 (6.0 
percent decrease) giving a modified value of P^ = 1.076. The value of 
LCCS decreases by approximately $315 or a relative change of 23 percent 
in the baseline value of $1369. 
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The information of Tables 6-1 through 6-5 can also be used to estimate 
the total uncertainty in life cycle cumulative savings due to uncertainty 
in different variables. If all the economic parameters are subject to 
variation a reasonable estimate of savings uncertainty can be obtained 
by the following: 


ALCCS 


prob 


N 


L j = 1 


/ 3LCCS 
\ '^j 


“’'j ) 


T i 
2 


(14) 


As an example, assume uncertainties of +10 percent in all fifteen of the 
variables listed in Table 6-1. The probable uncertainity estimate, using 
the data from the Table is: 


Toqus, Maine 
ALCCS 


prob 


= $1610 


This value is slightly larger than the present worth of cumulative savings 
of $1369 for Togus, given in Table 5.2-1. Had the probable uncertainity 
estimate greatly exceeded the cumulative savings, the risk of purchasing 
the solar system in anticipation of savings would have been greater, 
in direct proportion to the magnitude of the uncertainty in the indi- 
vidual variables. The results for the other sites are as follows: 

Albuquerque, New Mexico 
= $1941 

Cumulative Savings = $4536 
Forth Worth, Texas 

"■WSppPi, = $1114 

Cumulative Savings = $1096 
Madison. Wisconsin 

“'■“Sprob ' 

Cumulative Savings = $231 
Washington. DC 

'"■‘^“prob = *’906 

Cumulative Savings = $3038 
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TABLE 6-1 


UNCERTAINITY ANALYSIS FOR TOGUS, MAINE 


2 

Optimized Collector Area = 273 FT 


1 

i 

COST PARAMETER (x.) 

J 

NOMINAL 

VALUES 

NOMINAL 

VALUE 

DELTA 

3P1 

3P2 

’’‘j 

3LCCS 

3x. 

I 

ALCCS ! 

($) ! 
1 

(. 

AREA DEPENDENT COST (C^) 

12.750 

1.2750 

0.0 

0.0 

-313 

! 

-398 ! 

1 

AREA INDEPENDENT COST (C^) 

1413.000 

141.3000 

0.0 

0.0 

-1 

-162 I 

1 

FUEL COST (Cp) 

15.470 

1.5470 

0.0 

0.0 

457 

707 j 

DOWN PAYMENT/ INIT INV (D) 

0.200 

0.0200 

0.0 

-0.074 

361 

7 ^ 

FIRST YR. MI SC COST/ INIT INV (M) 

0.005 

0.0005 

0.0 

21.006 

-103090 

-52 ‘ 

FIRST YR. ASSESSED VAL/INIT INV (V) 

0.0 

' 0.0 

0.0 

0.0 

0 

0 

SALVAGE VAL/INIT INV (G) 

0.0 

0,0 

1 

0.0 

-0.196 

957 

0 

! ANNUAL MKT DISCOUNT RATE (d) 

0.085 

0.0085 

-286.35 

-7.987 

-37069 

-315 . 

I ANNUAL MKT RATE OF FUEL COST INC. (e) 

0.125 

0.0125 

\ 

252.55 

0.0 

67165 

1 

840 

• 

1 ANNUAL INT. RATE ON MORTGAGE (i) 

0.135 

I 0.0135 

0.0 

4.406 

-21563 

1 

-291 : 

i ANNUAL RATE OF GENERAL INFLATION (g) 
1 

0.100 

; 0.0100 

0.0 

0.954 

-4668 

-47 • 

j PROPERTY TAX RATE (t) 

0.0 

j 0.0 

0.0 

0.0 

° ! 

0 ! 

' EFFECTIVE INCOME TAX RATE (t) 

0.3 

1 0.0300 

0.0 

-0.838 

4099 

123 ' 

ANNUAL HEATING AND HOT WATER LOAD (L) 

22.800 

1 2.2800 
) 

0.0 

0.0 

310 

707 1 

ANNUAL SOLAR FRACTION (F) 

0.754 

j 0.0754 

0.0 

0.0 

9372 

707 


TABLE 6-2 
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C. UNCERTAINITY ANALYSIS FOR ALBUQUERQUE, NEW MEXICO 


Optimized Collector At^a = 168 FT^ 


j COST PARAMETER (x.) 

NOMINAL 

VALUES 

NOMINAL 

VALUE 

DELTA 

1 

3P1 

3x. 

3P2 

3x- 

J 

— 

3LCCS 

3x. 

J 

ALCCS 

($) 

AREA DEPENDENT COST (C^) 

12.750 

1.2750 

0.0 

0.0 

-192 

-245 I 

w 

; AREA INDEPENDENT COST (C^) 

1413.000 

141.3000 

0.0 

0.0 

-1 

-162 1 

1 FUEL COST (Cp) 

20.390 

2.0390 

0.0 

1 0.0 

425 

868 ‘ 

DOWN PAYMENT/ I NIT INV (D) 

0.200 

0.0200 

0.0 

-0.074 

262 

5 ; 

i FIRST YR. MISC COST/INIT INV (M) 

0.005 

0.0005 

0.0 

21.006 

-74889 

-37 

i FIRST YR. ASSESSED VAL/INIT INV (V) 

0.0 

0.0 

0.0 

0.0 ^ 

0 

0 ; 

i SALVAGE VAL/INIT INV (G) 

0.0 

0.0 

0.0 

-0.196 ■ 

695 

0 i 

ANNUAL MKT DISCOUNT RATE (d) 

0.085 

0.0085 

-285.35 

-7.987 I 

-65102 

1 

cn 

cn 

to 

ANNUAL MKT RATE OF FUEL COST INC. (e) 

0.125 

0.0125 

252.55 

0.0 1 

82459 

1031 1 

ANNUAL INT. RATE ON MORTGAGE (i) 

0.135 

0.0135 

0.0 

4.406 1 

t 

-15664 

-211 1 

ANNUAL RATE OF GENERAL INFLATION (g) 

0.100 

0.0100 

0.0 

0.954 \ 

-3391 

IJ 

-34 ; 

PROPERTY TAX RATE (t) 

0.0 

0.0 

0.0 

1 

0.0 [ 
1 

0 

0 1 

EFFECTIVE INCOME TAX RATE (t) 

0.300 

0.0300 

0.0 

-0.838 1 

2978 

89 , 

ANNUAL HEATING AND HOT WATER LOAD (L) 

17.200 

1.7200 

0.0 

< 

0.0 i 

504 

868 

! ANNUAL SOLAR FRACTION (F) 

0.931 

0.0931 

0.0 

0.0 

9318 

868 
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TABLE 6-3 


UNCERTAINITY ANALYSIS FOR FORT WORTH, TEXAS 


2 

Optimized Collector Area = 147 FT 


COST PARAMETER (x.) 

NOMINAL 

VALUES 

NOMINAL 

VALUE 

DELTA 

m 

8P2 

3x. 

0 

3LCCS 

3x. 

' 1 

ALCCS 1 

($) 

AREA DEPENDENT COST (C^) 

12.750 

1.2750 

0.0 

0.0 

-168 

-215 

AREA INDEPENDENT COST (C^) 

1413.000 

141.3000 

0.0 

0.0 

-1 

-162 

FUEL COST (Cp) 

13.010 

1.3010 

0.0 

0.0 

379 

493 

DOWN PAYMENT/INIT INV (D) 

0.200 

0.0200 I 0.0 

t 

0.0005 1 0.0 

1 

-0.074 

242 

5 

FIRST YR. MISC COST/INIT INV (M) 

1 

0.005 

21.066 

-69248 

-35 

FIRST YR. ASSESSED VAL/INIT INV (V) 

0.0 

0.0 

0.0 

0.0 

0 

0 

SALVAGE VAL/INIT INV (G) 

0.0 

0.0 

0.0 

-0.196 

643 

0 

ANNUAL MKT DISCOUNT RATE (d) 

0.085 

0.0085 

-286.35 

-7.987 

-26829 

-228 

ANNUAL MKT RATE OF FUEL COST INC. (e) 

0.125 

0.0125 

252.55 

0.0 

46817 

585 

ANNUAL INT. RATE ON MORTGAGE (i) 

0.135 

0.01 35 

0.0 

4.406 

1 -14485 

-196 

ANNUAL RATE OF GENERAL INFLATION (g) 

0.100 

0.0100 

0.0 

0.954 

i -3135 
! 

-31 

PROPERTY TAX RATE (t) 

0.0 

0.0 

0.0 

0.0 


0 

EFFECTIVE INCOME TAX RATE (t) 

0.300 

0.0300 

0.0 

-0.838 

i 2754 

83 

ANNUAL HEATING AND HOT WATER LOAD (L) 

19.360 

1.9360 

0.0 

0.0 

254 

493 

ANNUAL SOLAR FR.ACTION (F) 

0.736 

0.0736 

0.0 

0.0 

1 6692 

! 

493 



TABLE 6-4 


% 

UNCERTAINITY ANALYSIS FOR MADISON, WISCONSIN 


2 

Optimized Collector Area = 147 FT 


COST PARAMETER (x.) 

— 

NOMINAL 

VALUES 

NOMINAL 

VALUE 

DELTA 

3P1 

9X- 

J 

3P2 

3x. 

<J 

3LCCS 

ALCCS 

($) 

AREA DEPENDENT COST (C^) 

12.750 

1.2750 

0.0 


-168 

-215 

AREA INDEPENDENT COST (C^) 

1413.000 

141.3000 

.0.0 

|B 

-1 

-162 

FUEL COST (Cp) 

12.21 

1.221 

0.0 


332 

406 

DOWN PAYMENT/ IN IT INV (D) 

0.200 

0.0200 

0.0 

-0.074 

242 

5 

FIRST YR. MISC COST/INIT INV (M) 

0.005 

0.0005 

0.0 

21.066 

-69248 

-35 

FIRST YR. ASSESSED VAL/INIT INV (V) 

0.0 

0.0 

0.0 

0.0 

0 

0 

SALVAGE VAL/INIT INV (G) 

0.0 

0.0 

0.0 

-0.196 

643 

0 

ANNUAL MKT DISCOUNT RATE (d) 

0.085 

0.0085 

-286.35 

-7.987 

-17481 

-149 

ANNUAL MKT RATE OF FUEL CGST INC. (e) 

0.125 

0.0125 

252.55 

0.0 

38573 

482 

ANNUAL INT. RATE ON MORTGAGE (i) 

0.135 

0.0135 

0.0 

4.406 

-14485 

-196 

ANNUAL RATE OF GENERAL INFLATION (g) 

0.100 

0.0100 

0.0 

0.954 

-3135 

-31 

PROPERTY TAX RATE (t) 

0.0 

0.0 

0.0 

0.0 

0 

0 

EFFECTIVE INCOME TAX RATE (t) 

0.300 

0.0300 

0.0 

-0.838 

2754 

83 

ANNUAL HEATING AND HOT WATER LOAD (L) 

22.140 

2.2140 

0.0 

0.0 

183 

406 

ANNUAL SOLAR FRACTION (F) 

0.565 

0.0565 

0.0 

0.0 

7183 

406 










UNCERTAINITY ANALYSIS FOR WASHINGTON, DC 





flgurt 7-1. Economic Summary Chart for all Analyas Sitts 






7. SUMMARY AND CONCLUSIONS 

Solar energy is economically beneficial under the assumed economic con- 
ditions at Togus, Maine; Aibuquerque, New Mexico; Forth Worth i Texas; 
Madison, Wisconsion; and Washington, DC as shown in Figure 7-1. In 
Madison, Wisconsion, where the average solar insolation is 1191 Btu/ft^/ 
day and the conventional energy (electricity) cost is (0.042 $/kWh), a 
solar energy system of the type described in this report marginally 
profitable. Economic benefits from this solar energy system depend 
primarily on two factors: (1) maintaining or decreasing the initial 

investment required; (2) the continuing increase in the cost of convtfn- 
tional energy. The cap-'^ility to maintain or decrease the cost of the 
system relative to this present level is uncertain. It depends on favor- 
able tax treatment from the various levels of government, local through 
federal, as well as the continuing development of the solar energy in- 
dustry. On the other hand, increases in the cost of conventional energy 
are virtually assured. From the economic uncertainty analysis in Section 
6, where the conventional energy costs are medium to high, the savings 
with this system are 1.3 to 2.1 times more sensitive to increases in the 
conventional energy cost than to proportional increases in the solar 
energy system cost. This sensitivity serves to somewhat mitigate the 
risks. If the conventional energy costs are low, system cost increases 
and proportional increases in the cost of conventional energy equally 
impact the savings. 

The analysis and results given in this report can be used to guide a 
potential solar energy system buyer in evaluating the purchase of this 
type of DHW system. To do this, the solar insolation in the buyer's geogr- 
aphic area must be known. This data is available from several sources, 
including [5] and [6]. The cost of conventional energy must also be known. 
The local utility company can furnish rates from which a comparison cost 
based on 1000 kWh use can be computed in dollars per kWh. These values can 
then be compared with the characteristics of the analysis sites given in 
Section 3.1. The results for that analysis site can be ascertained from 
S'" iion 5.1 and 5.2. The primary economic parameters such as solar system 
costs, mortgage rates, inflation rates, discount rates, etc., are generally 
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known by the 'juyer for his area, bjviations in these econc/nic parameters 
from the values assumed in developing the results in this report can be 
evaluated from material included in Section 5. The aLCCS values given in 
Tables 6-1 through 6-5 were computed based fn a 10 percent increase in the 
economic parameter in question. A 10 percent decrease simply means changing 
the sign of the value in the appropriate table. Larger increases or decreases 
in an economic parameter cait also be obtained by multiplying the aLCCS value 
by the ratio of the desired increase to the 10 percent increase used in the 
original computation. 

As an example of the discussion above, assume the buyer has determined 
that the characteristics of his locale are similar to Fort Worth, Texas, 
and is considering the results reported for this solar energy system in 
Fort Worth. He notes that the reported savings from Table 5.2-1 is 
$1096, however, the conventional energy cost of his locale is $0.040/kWh, 
instead of the $0.044/kWh (Table 5.1-3) used in developing the Fort 
Worth saving. To modify the savings to consider the new rate the change 
is computed as: 


0»040 Q-P 4 4 ^ iQQ^ _ 9.1 ^(decrease) 

0.044 


In Table 6-3 for Fort Worth it can be seen that a 10 percent increase in 
fuel cost yields a value for aLCCS of $493. The impact on the Life Cycle 
Cost Savings of a 9.1 percent decrease in fuel cost can be computed as 
follows: 


aLCCS = * 492 = $449 

10.0 


Therefore, the new savings is: 

$1096 - $449 = $647 
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The buyer can evaluate the result of a change 1n any of the economic 
parameters In the same manner. However, he should be aware that the 
parameters are sometimes Inter-related and a change In one parameter 
may affect the aLCGS for several parameters. Consequently, the larger 
the changes the less the accuracy. However, approximate results may 
be obtained that prove of value In making a final decision. 


8 . 
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APPENDIX A 


F-CHART PROCEDURE 


APPENDIX A 


F-Cliart Procedure 

Modifications are made to f-Chart to enable the j,^*ogram to be used to 
perform economic analysis of the following: 

1. Systems that use heat pumps and fossil fuel space 
heating systems, as well as electric resistance 
heat. 

2 . Systems that use two different er<ergy sources fcsr 
•domestic hot water heating and space heating. 

The problem of analysis of the solar energy system with a conventional 
backup other than electric resistance heat is resolved by introducing 
a Coefficient of Performance (COP) whose value is dependent upon the 
type of backup system. Typical COP's of heat pumps are computed from 
a heat pump model which uses as inputs the ambient and building tempera- 
ture. Fossil fuel furnace COP's are assumed to be 0.60 unless different 
efficiencies, based on manufacturer's or other sources of data, are 
available. 

The problem of analysis with two different energy sources is resolved 
by adjusting the COP of the space heating system in accordance with the 
type of fuel used for the DHW system. This is necessary because the 
structure of f-Chart assumes electric energy to be the source for both 
space heating and domestic hot water. The adjustment factor is the 
adjusted ratio of the rates for the two energy sources used. The general 
expression for this is; 

CM _ DHW Auxiliary Fuel Rate ($/mi11ion Btu) 

"" SH Auxiliary Fuel Rate (i/million Btu) 

where the DHW Auxiliary Fuel Rate is the effective rate for 
fuel actually used and is equivalent to the electrical energy 
rate in a 100 percent efficient electrical hot water heater. 

The DHW auxiliary Fuel Rate will also be used for the value of 
Item Number 31 and 34 for systems of this configuration. 


SH COP 
or 

SH EfficiencyJ 


A-2 


The value of SH COP' is input to the modified f-Chart program. 

This value is used to compute an adjusted total load. The load, 
in turn, is used to derive the solar fraction which is input to 
the f-Chart economic analysis subroutine. 

Major considerations of the final report analysis procedure are the 
definitions of the loads that the system supports as it is analyzed 
in different geographic locations, and the sizing of the system to 
handle these loads at the various locations. The method it outlined 
in the following paragraphs. 

The monthly long-term heating load at the selected analysis sites is 
computed from the following equation: 

HL lt = UA*24*HDDlj - HTGEN 


where 


UA is the modified building energy loss coefficient 

HDD^j is the monthly long-term average heating degree days 

HTGEN is the internally generated heat computed from 
measured data. 

It is to be noted that UA is a modified parameter. The modification is 
to compensate for the fact that housing standards differ from U'^ation 
to location, i.e., the construction standards for a Florida house are not 
suitable for the New York environment. The UA factor used is derived from 
the ASHRAE 90-75 Standard [8] as a function of long-term heating degree 
days according to the appropriate U-value. The area, A, is derived from 
the building where the system is installed. 
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HTGEN is a factor that accounts for the part of the load whi({h is 
Internally generated. This is assumed to be the heat added which 
brings the building to the desired (comfortable) temperature when 
the outside ambient temperature is 65“F and no auxiliary heat is 
being added to the building. HTGEN, once derived, is assumed to 
be constant since it is a function of the life style of the occupants. 

The value of HL^j is the monthly long-term average heat load input 
to f-Chart. 

Additional technical and economic parameters that are input to f-Chart 
for the final report analysis are listed below with applicable 
comments . 

1. Air SH + WH = 1, Liq SH + WH = 2, Air or Liq WH Only = 3 

Comment; This is a definition of system type. The value 
is 1, if the system uses air collectors and supplies both 
space heat and domestic hot water; 2, if the system uses 
liquid collectors and supplies both space heat and domestic 
hot water; 3, if the system uses either type of collector 
and supplies only domestic hot water. 

2. (Flow rate/col. area) * (Spec, heat) 

Comment: If the system is an air system, this parameter is 

applicable. It is the air mass flow rate in Ib/min divided 
by the gross collector area multiplied by the specific heat 
of air at standard conditions. The value of this parameter 
is computed for the system at the actual installation site. 

This value is then maintained constant as the collector size 
is optimized for all analysis sites.* 

2 

*f-Chart uses an optimized value of 2.15 Btu/Hr °F Ft for this parameter 
In resizing a system, only the collector size is varied. The system is 
not given the benefit of further optimization. 
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3. 


eCmin/UA 


Comment; If the system is a liquid system and uses a liquid 
to air heat exchanger in the space heating loop, this parameter 
is applicable. It is the manufacturer's heat exchanger effec- 
tiveness multiplied by the minimum capacitance rate through 
the heat exchanger and divided by the building energy loss 
coefficient. If the heat exchanger effectiveness is unknown, 
a default value of 0.5 is specified. The capacitance, Cmin, 
is the minimum product of mass flow rate and specific heat, 
which usually occurs on the air side. The UA value 1s the mod- 
ified parameter applicable to the site. Deriving this value 
of UA has been previously discussed. The value of eCmin/UA 
is computed for the system at the actual installation site. 

This value is then maintained constant as the collector size 
is optimized for all analysis sites.* 

4. Collector Area 

Comment: This is the gross collector area which is optimized 
for all analysis sites. The optimization is extended to the 
actual installation site if an optimum sizing is not apparent 
in the original design. The predicted performance with optimal 
collector sizing is then compared to the predicted performance 
of the actual design and the actual measured performance. 

5. (ra) 

Comment: The basic value of (ta) was derived from the col- 
lector analysis program. This value is more consistent with 
actual operation than the manufacturer's or laboratory single 


*f-Chart uses an optimized value of 2.0 (dimensionless) for this parameter 
In resizing a system only the collector size is varied. 

The system is not given the benefit of further optimization. 
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panc.1 test values. If the system has a heat exchanger 
between collectors and storage, the derived value of 
(ta) was modified by the Fj^'/Fj^ factor as outlined 
in Section 2.4.4 of EES Report 49-3 (f-Chart Users 
Manual). [1] Note that the values input to f-Chart are 
assumed to be derived in accordance with ASHRAE specified 
method. 

6 . 

Comment: Same comment as Item 5. 

7. Incidence Angle Modifier 

Comment: In general, the default value of 0 is used. For 
evacuated tube collectors modeled as flat plate collectors 
the collector angle incidence modifier is obtained from the 
collector manufacturer. 

8. Number of Transparent Covers 

Comment: This is specified according to the characteristics 
of the collector. 

9. Collector Slope 

Comment: Collector Slope is changed according to the 
latitude of the site and the type of system. When the site 
analyzed is the existing site, the actual slope value is 
used. For other analysis sites the slope is computed as 
follows: 

• Latitude +10° if space heat and domestic hot water 

• Latitude if domestic hot water only 
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10. Azimuth Angle 


Conment; At sites other than the existing installation site 
the aximuth angle is 0°. At the existing site azimuth angle 
used for analysis was actual. However, any resulting per- 
formance degradation is noted. 

11. Storage Capacity 

Comment: This parameter is computed as the product of storage 
mass and specific heat divided by collector area for the exis- 
ting site. The same value of storage capacity is used for all 
sites. 

12. Effecting Building UA 

CoiTwent: The building UA, if not known, is derived from the 
measurement data contained in the Seasonal Report [7]. The 
computed value of UA is compared for reasonableness with a 
corresponding value of UA derived from ASHkAE Standard 90-75. 

For other analysis sites the value of UA is derived from ASHRAE 
Standard 90-75 as a function of building type and heating degree- 
days for each site. 

13. Constant Daily Building Heat Generation 

Comment; For residential type buildings, this parameter is 
derived from the measurement data contained in the Seasonal 
Report [7]. The derived value is held constant for all 
analysis sites. 

14. Hot Water Usage 

Comment: An effective average hot water consumption rate 
that accounts for actual load plus standby losses was 
computed from the following equation: 
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UUCP X UUAT 

HWCSMPEFF = Cp j TMAIN + TSETj * (TSET - TMAIN) * RHO J TMAIN + TSET ) 

Month 

15. Water Set Temperature 

Comment; The actual value of this parameter at the existing site 
is used for all analysis sites. 

16. Water Main Temperature 

Comment: The inputs for this parameter are a series of monthly 
values. The actual monthly value at the existing site is 
referenced to the average long-term ambient for the month for 
analysis at that site. For analysis at other sites the 
monthly value of TWAIN was established by site measurement 
at a nearby site referenced to the average long-term ambient 
for the month. (See Appendix C) 

17. City Call Number 

Comment: If the analysis site is located at a city listed in 

the November 1978 Input Data For Solar Systems that site is 
entered into the f-Chart data record. If the analysis site 
is not a part of the data record, an interpolative routine 
computes the data for any arbitrary site from nearby sites 
where data is available. 

18. Thermal Print Out by Month 

Comment: None 

19. Economic Analysis 

Comment: In general, all runs made for Final Reports specify 

print out of economic analysis. 
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20 . 


Use Optimized Collector Area - V, Specified Area 2 


Comment; In general the runs made for Final Reports use 
an optimized collector area. 

21. Solar System Thermal Performance Degradation 
Comment: A value of zero percent is used. 

22. -46. Economic Parameters 

Comment; The values of the economic parameter were worked 
out between MSFC and IBM for the Final Reports. The source 
of the value is given in the notes on page A-11. 


Residential 


Item 

Variable Description 

Value 

Units 

Source 

22 

Period of Economic Analysis 

20 

Yrs. 

SAI^ 

23 

Collector Area Dependent System Costs 



MSFC^ 

24 

Constant Solar Costs 



MSFC^ 

25 

Down Payment {% of Original Investment) 

20 

% 

SAI^ 

26 

Annual Interest Rate on Mortgage 

13.5 

% 

MSFC*^ 

27 

Term of Mortgage 

20 

Yrs. 

SAl' 

•1 

28 

Annual Nominal (Market) Discount Rate 

8.5 

% 

SAl' 

29 

Extra Insur. , Maint. in Year 1 

0.5 

% 

MSFC^ 


i% of Orig. Inv. ) 




30 

Annual % Increase in Above Expenses 

10.0 

% 

MSFC'^ 

31 

Present Cost of Solar Backup Fuel (BF) 



Actua r 

32 

BF Rise: %/Yr. = 1, Sequence of Values = 2 

1 
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Residential (Continued) 


Item 

Variable Description 

Value 

Units 

Source 

33 

Annual Rate of BF Rise 





Electricity 

12.5 

% 

MSFC'^ 

A 


on 

12.5 

% 

MSFC^ 


Natural Gas 

12.5 

% 

MSFC^ 

34 

Present Cost of Conventional Fuel (CF) 



Same as #31^ 

35 

CF Rise: tHr. = 1, Sequence of Values - 2 

1 



36 

Annual Rate of CF Rise 





Electricity 

12.5 

% 

MSFC'^ 

A 


Oil 

12.5 

% 

MSFC^ 


Natural Gas 

12.5 

% 

MSFC^ 

37 

Economic Print Out by Year = 1, 

2 


Analyst 


Cumulative = 2 



Option 

38 

Effective Federal State Income Tax Rate 





Residential 

30 

% 

SAl' 


Commercial 

48 

% 

MSFC^ 

39 

True Property Tax Rate Per $ of Original 

0 

% 

SAl' 


Investment 




40 

Annual % Increase In Property Tax Rate 

NA If m 

1 Is "O’ 

tl 

41 

Calc. Rt. of Return on Solar Investment? 



Analyst 


Yes = 1, No = 2 




42 

Resale Value (% of Original INvestment) 

0 


MSFC^’^ 

43 

Income Producing Building, Yes = 1, 



Site 


No = 2 



Dependent 

44 

Dprc.: Str. In. =1, Dc. Bal. = 2, 

2 

% 

MSFC'^ 


Sm-yr.-Dgt. = 3, None - 4 




45 

If 2, What % of Str. Ln. Dprc. Rt. Is Desired 

150 

% 

MSFC"^ 

rt 

46 

Useful Life for Deprec. Purposes 

20 

Yrs. 

MSFC"^ 


A-IO 


47. Economic COP for Auxiliary System 


Comment: This is a new parameter defined for f-Chart to 
account for economic analysis of solar systems having aux 
iliary backup other than electric resistance heat. The 
default values of this parameter are as follows: 

Heat Pump Auxiliary COP = 2 

Fossil Fuel Auxiliary COP = 0.6 

Electric Resistance COP = 1.0 


The value of the basic COP is modified, according to the method described 
on page A-2, to account for differences between the fuel used for the 
domestic hot water and the fuel used for space heating. 


A.-U 


NOTES; 


1. Source was Science Applications, Inc. (SAI) Draft Final Report 
on “Comparison of Solar Heat Pump Systems to Conventional 
Methods for Residential Heating, Cooling, and Water Heating," 
April 1979. 

E. These items were based on judgment and best experience. 

3. The actual current utility rates for the analysis sites 
selected were obtained. (See Appendix D). 

4. The assumption for final report analysis was that the backup 
system actually used for the installation was the same type 
of system that would be used If the solar system was not 
Installed. 

5. The declining balance technique never permits 100« depreciation 
of the asset no matter how long the period. The balance re- 
maining at the end of the system lifetime was treated, for 
accounting purposes, as salvage value. No other salvage 
value was presumed to exist. 


APPENDIX B 

ECONOMIC UNCERTAINTY ANALYSIS 
EQUATIONS 
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APPENDIX B 


ECONOMIC UNCERTAINTY ANALYSIS EQUATIONS 

1. Area dependent Investment costs (C^) 

ALCCSp/v = -PgA (%)(C^) 

2. Area Independent Investment costs (C^) 

ALCCS^,^ = -Pg (aC^)(C^) 

3„ Ratio of downpayment to initital investment (D) 


.LCCS„ = -(V . c,) ji - (-t) 


tf(N, 1, d) [' ■ f(N,'o, irlj 


4. Ratio first year's misc. costs to init. inv. (M) 


ALCCS 


M 


-(C^A + C^) [O - Ct) f(N» g, d)j (am) 


5. Ratio first year's assessed value to init. inv. (V) 


ALCCS, 


(c 


^A + C^) [t ( 1 - t) f(N. g. d)] (AV) 


6. Ratio salvage or resale value to init. inv. (G) 

r -1 


ALCCS 


G 


■(V + Cg) 


0 + d)^ 


(AG) 
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7. Annual market dlscoynt rate (d) 


ALCCS, 


8 . 


CpLF (1 - Ct) f(N, e, d) (Ad) 

"(V |f(N Jo?' lT 3d 

[^(1 - Ct) M + t (1 - t)v] 1^ f(N, g, d) - 

^ r'(N, 0, d) + 

- f^NToTiT) 3d ■ ji + d) 

^ 1^ f(N, 0, d) I (Ad) 

Annual market rate of fuel price increase (e) 


^ ff(N,^0, Ty 3d 


N+T 


ALCCS, 


CpLF (1 - Ct) 1^ f[h, e, d) (Ae) 


9. Annual Interest rate on mortgage (i) 


ALCCS. 


-(CjA + Cg) j(D - 1) (1 - t) : 

f(N, 0, d) 

flfTToTl) 

|rf(N, O; 1) - t (1 - D) [^1 - : 

1 

f(N, 0, 1) 

If f(N, 1. d) - t (1 - 0) f(a. 

i» d) 

' * f(N, 0. i) W "I 

1 (A1) 

li-3 



llitBMM 


10. Annual rate of general inflation (g) 


ALCCSg = -(C^A + Cg) [(1 - Ct) M + (1 - t) t 

Ig f (N, g, d) (Ag) 

11. Effective income tax rate (F) 


AUCS:^ = -CpLFCf{N, e, d) (At) 



(0-1 ) 


f^,_0, 

jar; o: 


+ (0 - 1) f(N, i. d) 


^ " fw; 0, i) 

- 


I f(N, 0, d) 


-t Vf(N, g, d) - C 

% 

(At) 


Mf(N, g, d) + 


12. Property tax rate (t) 

I 

ALCCS^ = -(C^A ■!- Cp) (1 - t) Vf(N. g, d) (At) 

13. Cost of conventional fuel in the first year (Cp) 


ALCCS(.p = P^LF (ACp)(Cp) 


14. Annual heating and hot water load (L) 
ALCCSp_ = P^CpF (aL) 


15. Annual load fraction supplied by solar (F) 


ALCCSp = P^CpL(AF) 


NOTE: 


Three functions used above require definition, as follows: 


f(N, a, b) 



|^f(N, a. b) = 


J 

b - a 




)( 


n + a 


I 


If f(N, a, b) 





APPENDIX C 

MONTHLY AVERAGE WATER 
SUPPLY TEMPERATURES 
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TABLE C-1 


MONTHLY AVERAGE WATER SUPPLY TEMPERATURES IN °F 


MONTH 

SITE NAME 

J 

F 

H 

A 

M 

J 

0 

A 

S 

0 

N 

D 

AVERAGE 

I 

T06US, ME 

43 

41 

42 

42 

52 

53 

62 

64 

62 

60 

55 

47 

52 ! 

ALBUQUERQUE, NM 

66 

66 

66 

70 

74 

76 

80 

83 

79 

74 

71 

66 

73 * 

1 

\ 

FORT WORTH, TX 

42 

49 

58 

65 

73 

80 

82 

83 

78 

63 

53 

49 

65 

1 

1 MADISON, WI - 

34 

37 

39 

50 

61 

68 

70 

72 

68 

63 

54 

36 

54 

f 

WASHINGTON, DC 

i 

42 

42 

52 

56 

63 

i7 

67 

78 

79 

68 

55 

46 

60 ; 




APPENDIX D 

ENERGY COSTS FOR 
ANALYSIS SITES 
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TOGUS, ME 

ELECTRICITY (RESIDENTIAL RATES) 


5.70 $/Month 
0.0446 $/kWh 

TAX 5% 

1000 kWh EFFECTIVE RATE = 0.053 $/kWh 


D~Z 


♦ O 


ALBUQUERQUE. NM 


r 


GAS 


(RESIDENTIAL) 


0-165 THERMS 0.0803/THERM 
165-340 THERMS 0.0826/THERM 
340+ THERMS 0.0966/THERM 

SERVICE CHARGE $1.25 EXAMPLE 

FUEL ADJUSTMENT $0. 21 14/THERM 30 THERMS * 0.2114 = 

TAX 4% 


ELECTRICITY 


(RESIDENTIAL) 


0-200 kWh 0.05294/ kWh 

200-800 kWh 0.04794/ kWh 

800+ kWh 0.03894/ kWh 

OR 

80C + kWh 0.04094/ kWh. 


1000 kWh EFFECTIVE 
NOV-MAY RATE = 0.069576 $/kWh 
YEAR-AROUND 

JUN-OCT 


FUEL RATE ADJUSTMENT $0.016680/ kWh 
SERVICE CHARGE $2.60 
TAX 4.5% 


FUEL OIL 


$0.999/G,AL+ 4% TAX 

I 


$6.34 


0-3 


FORT WORTH. TEXAS 


GAS 

0-1000 MCF 
1000-MCF 

SERVICE CHARGE 
TAX 

ELECTRICITY 

0- 25 kWh 
25+ kWh 
FUEL CHARGE 
SALES TAX 


$4.05/MCF 

$2.433/';CF 

0 

0 


$6.00 (MINIMUM) 
0.0285/ kWh 
$0.008899/ kWh 
4% 


1000 kWh EFFECTIVE RATE $0.0444/kWh 


FUEL OIL 


NOT USED IN FORT WORTH AREA 


MADISON. WI 


GAS 


0-2- 

20-50 

50+ 

THERMS 

THERMS 

THERMS 

$0. 28732/THERM 
0.27936/THERM 
0.26892/THERM 


ALSO 

FUEL RATE CHARGE 
TAX 

SERVICE CHARGE 

$0. 0762/THERM 

0. 

$2.00/M0NTH 

ELECTRICITY 

(RESIDENTIAL) 


i 1 

o o 
o 

100 kWh 
500 kWh 

$0.0360/ kWh 
0.0350/ kWh 



500-1000 kWh 0.0320/ kWh 
1000+ kWh 0.0275/ kWh 

ALSO FUEL RATE CHARGE (JAN) $0.00607/kWh 
TAX 0 

SERVICE CHARGE $2.00/M0NTH 
1000 kWh EFFECTIVE RATE $0.04167/kWh 
FUEL OIL 


$0. 919/GAL 

TAX 0 FOR RESIDENTIAL 4% FOR COMMERCIAL 
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WASHINGTON, DC 


GAS 


$5.00/M0 SERVICE CHARGE 1 THERM » 100,000 Btu 

$0. 3255/THERM + 5 % TAX 


ELECTRICITY 

(RESIDENTIAL RATES) 





$5.00/M0 SERVICE 

CHARGE 






NOV - 


« 


• . 

JUNE - 

■ OCT 


WINTER RATES 





SUMMER 

RATES 


0 - 600 

kWh 

0.06024 

$/ kWh 


0 - 600 

0.06024 

$/kWh 

600 - 1500 

kWh 

0.05334 

$/kWh 


600 - 1500 

0.06924 

$/ kWh 

1500 + 

kWh 

0.04289 

$/kWh 


1500 + 

0.26638 

$/ kWh 


TAX 16'/ OF FIRST $15.00 ($2.40 MAX) 

FUEL CHARGE 0.01500 $/ kWh (INCLUDED IN ABOVE RATES) 

1000 kWh EFFECTIVE RATE = 0.0G75 $/ kWh YEAR-ROUND 

FUEL OIL 

$0.989/6AL + TAX 5'; 
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